Acceptance corrections to net baryon and net charge cumulants 
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We show that the effect of finite acceptance drastically influences the net-baryon and net-charge 
cumulants, which are believed to be sensitive probes of the QCD phase diagram. We derive the 
general formulae that relate the true cumulants K n which reflect the full dynamics of the system 
with the actually measured cumulants c„ for a given acceptance, modeled by a binomial probability 
parameter p. We find that this relation involves additional moments which cannot be expressed 
by cumulants and should be measured in order to extract any potential information about the 
QCD critical point. We demonstrate that for a wide range of the true cumulant ratios K n /K m the 
measured ratios c„/c m quickly converge if p < 1/2, which makes the interpretation of the data very 
challenging, especially in case of the net-proton cumulants. Our study further suggests that the 
measurement of net-charge cumulants may be more advantageous for the investigation of the QCD 
phase diagram. 

PACS numbers: 25.75.-q, 24.85.+p, 24.60.-k 



I. INTRODUCTION 



The exploration of the phase structure of Quantum Chromodynamics (QCD) has been a central topic in the study 
of the strong interactions for many years. Theoretically, the QCD phase diagram has been investigated with either 
effective models or by first principle Lattice QCD calculations. Experimentally, the phase structure of strongly 
interacting matter is studied with the help of relativistic heavy ion collisions. Experiments at the Relativistic Heavy 
Ion Collider (RHIC) and recent results from the Large Hadron Collider (LHC) have indicated that the collisions at 
the very highest energies produce a strongly interacting quark-gluon plasma at nearly vanishing net-baryon density. 

Meanwhile Lattice QCD (LQCD) calculations have established that for vanishing net-baryon density the transition 
between the hadronic and partonic "phase" is an analytic cross over [l| with a pseudo-critical temperature of T c ~ 
160 MeV [2|,l3j|- Therefore, true phase transitions, if at all present, will be located at finite values of the baryon density. 
Indeed, model calculations predict a first order coexistence region at high baryon density (several times nuclear matter 
density) and moderate temperatures, which ends in a critical point (for a compilation of various model results see, e.g., 
Q). The precise location of this critical point, however, is not yet known from theory. Model calculations typically 
predict the location of the critical point to be at rather large values of the baryon number chemical potential, fiB, 
I which currently cannot be explored by Lattice QCD methods, due to the fermion sign problem. 
• • ' Experimentally, the region of the phase diagram at finite density and moderate temperature is explored by heavy 
ion collisions at various energies. Freeze out points extracted from experiments range from values of nearly vanishing 
baryon chemical potential for LHC and top RHIC energies to values of hb — 500 MeV at y/s ~ 5GeV [5j. Since 
the location of the critical end point (CEP) is not really known the entire accessible region of the phase diagram 
needs to be explored by scanning the full range of available beam energies. Both the CEP as well as the first-order 
phase transition are associated with characteristic fluctuations - long range for the second-order transition at the 
CEP, and possible spinodal instabilities in case of a first order transition 0, 0|. Therefore fluctuations of various 
quantities such as particle multiplicities or mean transverse momentum have been proposed as suitable observables 
Q. Measurements of this sort have been carried out by the NA49 collaboration [9|, [lO| at the CERN SPS and are 
also part of the recently started RHIC beam energy scan The first results by the NA49 collaboration |jj licj . 

which concentrated on the variances of various particle ratios and the transverse momentum, showed only very small 
deviations from the expected Poisson fluctuations of a Hadron Resonance Gas (HRG). 

Meanwhile it has been realized that higher order cumulants would be more sensitive to the fluctuations associated 
with the second order transition, including the CEP [12l4l4| . The principle reason is that they scale with a higher power 
of the correlation length, which will be finite due to the rather short lifetime of the system created in these collisions. 
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While the cumulants of many distributions may carry information about the CEP, theoretically the cumulants of 
the net-baryon number and of the net-charge are preferable since they are well defined as derivatives of the QCD 
partition function with respect to the appropriate chemical potentials. In addition, fluctuations of conserved quantities 
are less affected by final state interaction in the hadronic phase [l5j . Of course in an actual experiment, where the 
baryon number and electric charge of the entire system are conserved, corrections need to be applied to be able to 
compare with theoretical calculations which are typically carried out in the grand canonical ensemble, where charges 
are conserved only on the average. This point has been recently addressed, e.g., in [llij . 

In addition to corrections due to charge conservation, in a real experiment one only observes a fraction of the 
final state particles. In case of the net-baryon distribution, typically the neutrons are not observed. Furthermore 
there are acceptance cuts as well as efficiency corrections for the protons, so that less than 50% of baryons are 
actually observed. As we shall argue in this paper, this will lead to substantial corrections for the observed cumulants 
that need to be taken into account before any conclusions about critical fluctuations can be drawn. Specifically we 
will show that it is not sufficient to know the distribution of net-baryons (or net-charges) P{Nb — Ng) in order to 
predict the observed cumulants of the net-proton distribution. Instead, information for the full baryon/anti-baryon 
distribution, P(Nb, N§), is needed, which is not easily calculated in theory 1 . The same is true for the net-charge 
distribution, however in this case corrections are smaller since a larger fraction of charge particles are observed in a 
typical experiment. On the positive side, we will show that by measuring various factorial moments of the observed 
distribution one is able to extract the cumulants of the physical net-baryon distribution, P(Nb — Nq), which is, after 
all, the desired goal. 

We note that the question of acceptance corrections to various moments of particle distributions have already been 
discussed in the literature. For example corrections to the variance or second order cumulants of the net charge 
and net baryon number distributions have been investigated in Refs. [13413, while acceptance effects on factorial 
moments have been addressed in [20Ll2l|. Here we will concentrate on the corrections to higher order cumulants which 
are central to the discussion of the QCD critical point. 

This paper is organized as follows. In the following section we will derive the relation between the cumulants of the 
underlying distribution to that of the actually measured distribution, and we will show that additional information 
other than cumulants is required. In the next section we will present some examples to illustrate the problems which 
arise if the acceptance corrections are large, as it is the case for the measurement of the net-baryon distribution. 
Before we conclude we will discuss a few observations pertinent to our results. 



II. CUMULANTS 



Suppose we have an underlying probability distribution P(JVi, jVa) which captures the full dynamics of the system. 
Subsequently, we will refer to P(Ni,N2) as the distribution within the "required acceptance", where by required 
acceptance we mean that all particles that are necessary to capture the relevant physics are measured. This does not 
necessarily imply that all particles in the final state need to be measured. 2 Here Ni and N2 stand for either baryons 
and anti-baryons, i.e., N\ = Nb and N2 = N B or positively and negatively charged particles, i.e., iVi = N + and 
N2 = N-. Let further K n be the net-baryon or net-charge cumulants associated with the underlying distribution 
P(Ni, N2). Of course in a real experiment additional cuts such as in the transverse momentum or rapidity will have to 
be applied. In addition, the efficiency to detect particles will be smaller that 100% (in case of neutrons it is typically 
zero). To reflect this let us introduce the probability distribution p(ni, 112) of the actually observed multiplicities m 
and H2- Let us further denote the net-baryon or net-charge cumulants of p{n\ 1 n2) by c n . Finally, we assume that all 
acceptance corrections may be modeled by a binomial probability distribution 



00 00 ( 

p(m,n 2 )= J2 E P (Ni,N 2 ) rr pFCl-pi) 
^— ' ' ni (iVi — n\ ! 

JVi=rai N 2 =n 2 V ; 
WAT \\ P 2 I 1 P2j 

n 2 !(A'2 - n 2 y. 



1 In principle P(Ng , iVg) can only be calculated for systems where a quasi-particle description applies. 

2 To which extent an actual experiment does have the "required acceptance" is not easy to tell and naturally depends on the specific 
physics under consideration. 

3 In reality, more sophisticated models than the binomial distribution are likely required to properly address the various acceptance 
corrections (see e.g. Il8l0 . However, using the binomial distribution will illustrate the essential issues related with acceptance corrections. 
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The parameters p\ and p 2 describe all possible acceptance effects in our system. For example, the detector accep- 
tance is naturally modeled by the binomial distribution and it can be different for different particles. Additional cuts 
in the transverse or/and longitudinal momenta also introduce acceptance correction. In the context of the cumulants 
of the net-baryon distribution, it has been argued in Refs. [22I |23 | that the fact that neutrons are not observed in 
a typical experiment may also be modeled by a binomial distribution. Thus, Eq. ([T]) applies as well. Finally, we 
note that the different sources of limited acceptance can be represented by an effective parameter, which is simply 
the product of all binomial probabilities. For example, in case of the net-proton distribution: n\ = n p (number of 
protons), n 2 = rip (number of anti-protons), see Eq. ([1]), and 

Pl=p?-p?-pY..., (2) 

where pf represents the fact that only protons are measured instead of baryons (pf ss 1/2), pY characterizes the 
cut in momentum for measured protons, p^ is a detector efficiency for protons in a given momentum cut, and so on. 
Similar relation holds for p 2 . 

In order to proceed it is convenient to introduce the factorial moments: 

/N^. N 2 \ \ ^ ~ N \ N l 

Fa s ( m - o. w - tn ) - g, E m ' ■ V. - o- m - w- ■ (3) 

fik=(j ^ttt ^tt ) = V V p(n 1 ,n 2 ) 1 ^7 '-j-rj. (4) 



Using Eq. ([T]) it is straightforward to show (see e.g. [20|, |2l|): 

hk=p\-P2- F lk - (5) 

The next step is to express the cumulants K n of the distribution P(N% — N%) in terms of the factorial moments Fn~. 
Using the relation, Eq. ((5]) , we can then write K n in terms to the measured factorial moments fn~ for given acceptance 
parameters p\ and pi- Finally can be expressed by the cumulants c„, which will allow us to relate K n to c m . 
However, as we shall see below, the factorial moments fa- and Fik in general cannot be solely expressed in terms of 
the cumulants of the net distributions, c„ and K m , respectively. Additional terms, not related to cumulants will arise. 
It is useful to define the moment generating function: 



h(z) = V p(,S) 



z 5 



8 1 

00 00 

= E E P(Ni,N 2 )[l-(l-z) Pl ] Nl [i-(i- z -i) P2 ] N \ ( 6 ) 

N 1 =0N 2 =Q 

where 5 = rii — n 2 [see Eq. (fTJ)] and p{5) is the net multiplicity distribution 

P( S ) = E P( n l> n 2) S ni-n2-8- (7) 

71 1 ,712 

Equation §§§ allows to calculate the cumulant generating function 

00 £ 

9 (i)=ln[/ l (e t )]=^c fe ^. (8) 
fe=i 

For P i = p 2 = 1 in Eqs. (|6l8p we obtain the moment- and the cumulant generating function, G(t) = g(t)\ pi=P2=1 , of 
the net multiplicity distribution P(Ni — N 2 ). By definition, the cumulants c„ and K n read: 

d n g(t) d n G{t) 

The above equations allow us to calculate the cumulants c n and K n . To obtain the relation between them is 
straightforward but tedious. For instance both c 2 and K 2 contain Fn [see Eq. ©], so we express Fu by K 2 and 
substitute into the expression for c 2 . In this way we can relate cumulants c„ by cumulants K m and the factorial 
moments Fik- 
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First, let us relate the measured cumulants c„ with the cumulants K m of the underlying distribution P (N\ — N 2 ). 
Here we present the results up to C4 and we assume p = P\ = p 2 - 111 the Appendix we give the general relations up 
to eg [HI (see also [25l. [2^|) with an arbitrary p\ and p 2 . We obtain 

ci= P Kx, (10) 

c 2 = p(l-p)N+p 2 K 2 , (11) 

c 3 = p(l - p 2 )K x + 3p 2 (l - p) (F 20 - F Q2 - NK X ) + p 3 K 3 , (12) 

and more complicated: 

c 4 = Np(l - p) - 3iVV(l - P? + 6p 2 (l - p){F m + F20) - 12AV(1 - p)(A 20 - F 02 ) 
+6Ap 3 (l - p){Kl - K 2 ) + p 2 (l - p 2 ){K 2 - 3A 2 ) 

+6p 3 (l - p)(Fo3 - F 12 + F Q2 + F 20 - F 2 i + F 3Q ) +p 4 Ki. (13) 

To simplify the notation we have introduced 

N= (Ax) + (N 2 ) =F 10 +F 01 . (14) 

As can be seen from the above equations it is impossible to relate the cumulants c n solely with the cumulants K m 
and the parameter p. In Eq. (fTTj) the value of N is present. In c 3 the information about F 2 q = (N\(Ni — 1)) and 
Fq 2 is needed. In Eq. (fT3")) A30, F03 and mixed factorial moments F12 and F 2 i appear. As already discussed in 
the introduction, these additional terms would have to be calculated in a given theory, in order to make a reliable 
predictions for the measured cumulants c„. In general this is difficult and likely can only be done if a quasi-particle 
description for the dynamics applies. 

It is interesting to consider two limits of p: for p — 1 we obviously obtain c n = K ni and for very small p the 
following relations hold 

C2 ~ (ni) + (n 2 ) , c 3 « ci, c 4 ~ c 2 , (15) 

where (m) — p(Ni). In the limit p — > the ratios of cumulants read: C3/C1 = 1, c^jc 2 = 1. It means that even if 
the cumulant ratio K n /K m contains information about the QCD phase diagram (or any other interesting physics) 
this information is lost if the acceptance is too small. In other words in the limit p^O the Skellam distribution is 
obtained, i.e., we are in the limit of Poisson statistics. 

To summarize, Eqs. (|10H13[) allow to obtain the values of cumulants c n at a given acceptance parameter p, knowing 
the cumulants K n and certain factorial moments Fik of the underlying distribution P (Ni^N 2 ). This distribution 
according to our assumptions is obtained for the "required acceptance", defined by p = 1. 

In practice, the more interesting question is to which extent one can express the cumulants K n , which contain the 
relevant physics, in terms of the measurable cumulants, c m , and factorial moments, fik- Actually this can be easily 
done using Eqs. (I10H13[) and the relation ([5]): 

pKi = a, (16) 
p 2 K 2 = c 2 -n(l-p), (17) 
p 3 K 3 = c 3 - ci(l - p 2 ) - 3(1 -p){f 20 - fo2 -net), (18) 

and 

p 4 K 4 = c 4 - np 2 (l - p) - 3n 2 (l - p) 2 - 6p(l - p)(/ 20 + fa) + 12 Cl (l - p)(fa - fa) 
-(l- P 2 )(c 2 - 3c?) - 6n(l - p){c\ ~ ca) 

-6(1 -p)(/ 03 - /12 + fa + /20 - hi + fa)- (19) 
where we introduced the notation 

n=<m) + <na) = /io + /oi. (20) 

These equations can be directly used to extract the values of the true cumulants K n that characterize the system 
with p = 1, which we assumed to capture the relevant physics. We repeat: The measurement of the cumulants c m is 
not sufficient to extract the desired cumulants K n . The additional measurement of various factorial moments fa is 
required as well. 

For example, in the context of the net-baryon fluctuations, K n denote the net-baryon cumulants. To determine 
K n all we need as input is the value of the acceptance parameter p, the measured cumulants c„ of the net-proton 
distribution, and the factorial moments fa measured for protons and anti-protons. 
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III. EXAMPLES 

In the following we present a few examples where we show the relation between the observed cumulants c„ and 
the cumulants of the underlying distribution, K n . These are merely examples to illustrate the situation and certain 
assumptions about the necessary factorial moments will have to be made. 

A. C3/C1 

As an example we show the dependence of the ratio C3 jc\ as a function of the acceptance parameter p for different 
values of K3/K1. As seen from Eq. (fT2]) . C3 depends only on cumulants K\ and A3 and factorial moments F20 and 
F 02 - To proceed we denote 

(iV?>-<^> 2 = W(l + a): (21) 

where i = 1,2 and a is a free parameter that allows to change the widths of the distributions P(Ni) and P (N2) 4 - 
Substituting Eq. fl5T) to Eq. JEJ) and using Eq. ([TOjl. we obtain: 

^ = l-p* + 3ap(l-p)+p 2 ^. (22) 
ci Ai 

In the context of the net-baryon fluctuations it is expected that the width of the baryon (anti-baryon) distribution 
is comparable to the Poisson distribution [ll|, that is \a\ << 1. However, when net-charge fluctuations are studied 
the parameter a is expected to be of the order of one 5 . In Fig. Q]and Fig. [5] we present the dependence of c^/ci as a 
function of p for four values of K3/K1 = —1, 0, 0.5 and 1, and we note that, as previously discussed, A3/A1 equals 
to the value of C3/C1 at p = 1. In Fig. [I] we assume that the multiplicity distribution is narrower than the Poisson 
distribution: a = —0.1 in the left plot, and a — —0.5 in the right plot. 




Figure 1: The measured cumulant ratio c^/ci as a function of the acceptance parameter p for four values of K3/K1 = —1, 0, 
0.5 and 1. K3/K1 equals czjc\ at p = 1. In the left plot a = —0.1 and in the right plot a — —0.5, see Eq. (|21[) . For the 
measurement of the net-proton cumulants at STAR the realistic value of p is smaller than 1/2 and arguably close to 1/5. 

In Fig. [5] we choose a to be positive, i.e., the multiplicity distribution is broader than the Poisson distribution: 
a = 0.5 in the left plot, and a = 1 in the right plot. 

In case of net-baryon fluctuations the fact that we only measure protons (anti-protons) instead of all baryons (anti- 
baryons) already introduces p < 1/2. In addition finite detection efficiencies and cuts in the transverse momentum 
reduce the value of p even more. We estimate that for the STAR measurement [llj the parameter p m 1/5 may not be 
far from reality (assuming that the cut in rapidity is good enough to capture the relevant physics, and its contribution 
to p can be neglected). As seen in Figs. [1] and [2] already the value of p = 1/2 shifts the observed ratio C3/C1 into 



4 For simplicity we assume that a is the same for Ni and N2 multiplicity distributions. 

5 The charge multiplicity distribution in heavy-ion and proton-proton collisions can be described by the negative binomial distribution, 
see, e.g., |27|| . 
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Figure 2: The same as in Fig. [T] except in the left plot a = 0.5 and in the right plot a = 1, see Eq. (|21[) . 

positive values even if K3/K1 is negative. Moreover, for different values of K3/K1 the values of c^/ci quickly converge 
which makes the interpretation of the experimental results challenging. Taking a more realistic value of p = 1/5 we 
see that it is practically impossible to distinguish between different values of K3/K1, unless K3/K1 << — 1. 

B. c 4 /c 2 

It is interesting to repeat the same exercise for C4/C2. However, in this case we need additional information about 
-F3O) F03 and the mixed factorial moments F21 and F\ 2 . We investigated different values of parameters and we obtain 
similar plots as in Figs. [T]and[5J To present a few examples, we will assume that i 7 ^ = Fm 6 , and introduce the 
following parametrization 

F 2Q = (Nif + aiNx), (23) 
^30 = (Nif + 3a (Nif + 2a 2 (Ni) , (24) 
F21 = F 20 (N^ + (3 (N^ 2 . (25) 

Equations (|23|24| result from the assumption that P(N\) is given by the negative binomial distribution. The parameter 
a is identical as in Eq. (|2ip , and controls the width of the multiplicity distribution. f3 reflects the correlation between 
N\ and N%. To further reduce the number of parameters we assume that 7 : 

K 2 = 2 (Nf) - 2 (N X N 2 ) w 2 (A^ 2 ) _ 2 (A^) 2 . (26) 

Finally, the value of is determined from the ratio K4/K2 that is an input in our calculation. 

In Fig. [3] and Fig. [3] we present the dependence of C4/C2 as a function of the acceptance parameter p for five values 
of K 4 /K 2 = -5, -1, 0, 1, and 5. In Fig. Q]we assume N = 100, see Eq. jl4]), and a = -0.1. In the left plot (3 = 0.01, 
and in the right one f3 = —0.01. 

In Fig. 0]we assume N = 100 and a = 0.5. In this case the multiplicity distribution is broader than Poisson, which 
is expected if the net-charge cumulants are investigated. In the left plot (3 = 0.01, and in the right one j3 — —0.01. 

As seen from Figs. [3] and 0] the convergence of C4/C2 for various values of K4/ K 2 is even more rapid than in case 
of C3/C1. Taking p — 1/5, as in the STAR measurement of net-proton cumulants, we see that even if there is any 
information about the QCD phase diagram in the net-baryon cumulants, this information is strongly diluted if only 
protons and anti-protons are measured. 8 At p = 1/5 all lines are practically indistinguishable, unless \K±/K.2\ » 5 
(it should be of the order of 50). 



We note that this choice corresponds to a system at zero baryon chemical potential in the context of the net-baryon cumulants. However, 
the formalism presented in this paper is not restricted to this choice. Also, the corrections shown here are expected to be similar in case 
of finite chemical potential. 

We checked that modifications of this assumption do not introduce any new qualitative features. 

Our finding naturally explains the large difference between net-baryon and net-proton cumulants observed in the UrQMD model \2>\\ . 
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Figure 3: The measured cumulant ratio C4/C2 as a function of the acceptance parameter p for five values of K4/K2 ~ — 5, — 1, 
0, 1, and 5. K4/K2 equals C4/C2 at p — 1. In both plots a = —0.1. In the left plot (3 = 0.01, and in the right one f3 = —0.01, 
see Eqs. (|23II25[) . For the measurement of the net-proton cumulants at STAR the realistic value of p is smaller than 1/2 and 
arguably close to 1/5. 




Figure 4: The same as in Fig. [3j except in the left plot a = 0.5, /3 = 0.01, and in the right plot a — 0.5, /3 = —0.01. 

IV. COMMENTS 

Let us discuss a few observations pertinent to the results from the previous sections: 

(i) We demonstrated that a limited acceptance 9 makes a direct interpretation of measurements of higher order 
cumulants very difficult. Instead, we propose that the measurement of additional factorial moments, as shown in 
Eqs. (|16H19p . or the more general ones in the Appendix, is the way to proceed. This allows a rather unambiguous 
determination of the true cumulants K n . 

(ii) The main problem with the measurement of the net-proton fluctuations is that the maximal value of p is 
approximately 1/2. 10 Moreover, additional cuts in the transverse momentum are usually needed. Also detection 
efficiencies are never 100%. As a result, p < 1/2 which makes the interpretation of the data very challenging. Even 
if one measured the additional factorial moments as proposed in this paper, their contribution becomes large for 
small values of p [see Eqs. ()16til9|) ] and thus need to be determined very precisely. This clearly favors the analysis 
of net-charge fluctuations, where there is no problem equivalent to the neutrons. Thus the parameter p will likely be 
larger than 1/2 and the extrapolation to the full acceptance by measuring the necessary factorial moments will be 
much more reliable. The cumulants of the net-charge fluctuations should also be sensitive to the CEP, however with 



In particular the difficulty in measuring neutrons in the context of net-baryon fluctuations. 

Of course if one compares with theoretical predictions which directly calculate the proton cumulants, such as [2S(. the effective binomial 
probability p is closer to one, similar to that for charged particles. However, it is not clear how to determine the proton cumulants in a 
model independent way. 



a potentially smaller overall magnitude [l2j ■ 

(iii) For the purpose of this paper, we have assumed that the "required acceptance", i.e., p = 1, captures the 
relevant physics. To which extent this is the case for an actual experiment, such as STAR, is not clear and difficult 
to estimate. For example, it is not clear to us what the size of the rapidity window should be, in order to capture 
the relevant physics. One estimator would be the width of a thermal fireball, which is of the order of one unit of 
rapidity. However, it may very well be larger, and, therefore, the true value of p may very well be smaller than 1/5 
in the STAR measurement. 

(iv) Given the expression for the true cumulants K n in terms of the measured cumulants c m and factorial moments, 
Eqs. (|16H19p . it may be conceivable that one could explore the necessary rapidity range experimentally. As one 
increases the rapidity window the resulting K n should approach an asymptotic value once the window captures the 
relevant physics. On the other hand, doing the same exercise with the measured cumulants c n only may lead to false 
conclusions, especially in case of the net-proton distribution. Since the binomial parameter p is already very small 
reducing the rapidity window will only lead to small variations (see Figs. [T]|4]) suggesting that the physics does not 
change. 

(v) In our examples for C4/C2 we find that the choice for the correlation term ~ j3 in Eq. (|25|) leads to significant 
corrections for the measured cumulant ratio. Therefore, these correlation terms need to be measured precisely and 
we expect them to play an even stronger role in higher order cumulant ratios, such as cg/c2- 

(vi) We note the formalism presented here may be applied to higher order cumulants, such as Kq, in a straightforward 
but tedious way (see Appendix). We believe that the problem to resolve the true cumulants gets even more difficult, 
the higher the order of the cumulant under consideration is. 

(vii) Instead of measuring the cumulants one may consider other moments which are independent of the (binomial) 

acceptance corrections. For example the scaled factorial moments (l8l. [20L l2l| . — ^- rfe p = ^ ns , see Eqs. 

QUI), would be such an alternative. Or one could construct combinations of moments similar to the one suggested 
in [19( for the second moments. While these alternative moments should in principle carry similar information as 
the cumulants, they are not fully determined by the net-baryon distribution, P(Nb — Ng). Therefore, a theoretical 
determination of these moments, e.g. from Lattice QCD, will necessarily require some mo del- dependent assumptions. 

V. CONCLUSIONS 

In this paper we have studied the effect of finite acceptance on higher order cumulants of net-baryon and net-charge 
distributions. To this end we have folded the true probability distribution with binomial distribution in order to 
simulate the finite acceptance corrections. Our main finding is that the task of extracting the cumulants K n of the 
true distribution requires the measurement of not only the cumulants c n of the measurable distribution but in addition 
of various factorial moments, which cannot be expressed in terms of cumulants. We also demonstrated that for various 
values of the true cumulants ratios K n /K m the measured ratios c n /c m quickly converge with decreasing acceptance 
parameter p. This makes the physical interpretation of the net-proton cumulants very difficult. We further argued 
that it may be advantageous to investigate the cumulants and factorial moments of the net-charge distribution, since 
in this case acceptance corrections are considerably smaller than for the net-baryon distribution. This may allow for 
a more reliable extraction of the true cumulants via the methods presented in this paper. 
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Appendix A: General relations 

Here we present the general relations between cumulants K n characterized by the "required acceptance", and the 
measurable factorial moments at a given acceptance parameters p\ and p2- As seen in Eqs. (|17H19j) it is not possible 
to express cumulants K n solely by cumulants c m but also factorial moments fa appear. Therefore, here we express 
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K n solely by the factorial moments fik- This significantly simplifies our notation but makes no difference for an 
experimental application. We obtain 

K x = (Aq) - (N 2 ) , (Al) 
K 2 = N-K 2 + F 02 -2Fu + F 2a , (A2) 
K 3 = K\ + 2K l — F03 — 3Fo2 + 3iq 2 + 3F 2 q — 3F 2 i + F 3 q 

-3K 1 (N + F 02 -2F 11 +F 20 ), (A3) 
K A = N- 6Kf + F 04 + 6F 03 + 7F 02 - 2F n - 6F 12 - 4F 13 

+ 7.T20 — &F 2 i + 6F 22 + 6-F30 — 4i^3i + i*4o 

+ \2Kl{N + F 02 - 2F n + F 20 ) - 3(N + F 02 - 2F n + F 20 ) 2 

- 4Ki(K\ - F 03 - 3F 02 + 3F 12 + 3F 20 - 3F 21 + F 30 ), (A4) 

and K5 and Kq are more complicated 

K h = K l + 2AKf - F 05 - 10F 04 - 25F 03 - 15F 02 + 15F 12 + 20F 13 + 5F U 
+ 1 5F 20 - 15F 21 - 10F 23 + 25F 30 - 20F 31 + 10F 32 + 10F 40 - 5F 41 + F 50 

- 60Kf(N + F 02 - 2Fu + F 20 ) + 301^ (N + F 02 - 2F n + F 20 ) 2 
+ 20Kf(K 1 - F 03 - 3F 02 + 3F 12 + 3F 20 - 3F 21 + F 30 ) 

- 10(N + F Q2 - 2F n + F 20 )(ifi - F 03 - 3F 02 + 3F 12 + 3F 20 - 3F 21 + F 30 ) 

- 5K!(N + F 0i + 6F 03 + 7F 02 - 2Fu - 6F 12 - 4iq 3 + 7F 2a - 6F 21 + 6F 22 

+ 6F 30 -4F 3 i +F 40 ), (A5) 

and 

K 6 = N- UOKf + F 06 + 1 5F 05 + 65F 04 + 90F 03 + 31 F Q2 - 2F n - 30F 12 - 80F 13 
- 45F 14 - 6F15 + 31F 20 - 30F 21 + 30F 22 + 30F 23 + 15F 24 + 90F 30 - 80F 31 
+ 30F 32 - 20F 33 + 65F 40 - 45F 41 + 15F 42 + 15F 50 - 6F 51 + F 60 
+ 360Kf(N + F 02 - 2F n + F 20 ) - 270K 2 (N + F 02 - 2F n + F 20 ) 2 



+ 30{N + F 02 - 2Fu + F 20 ) 6 - 120Kf{K 1 - F 03 - 3F 02 + 3F 12 + 3F : 



20 

- 3F 21 + F 30 ) + 120Xi (N + F 02 - 2Fu + F 20 )(ifi - F 03 - 3F 02 + 3F 12 

+ 3F 20 - 3F 21 + F 30 ) - 10(^i - F 03 - 3F 02 + 3F 12 + 3F 20 - 3F 21 + F 30 ) 2 
+ 30Kf{N + F 04 + 6F 03 + 7F 02 - 2F n - 6F 12 - 4F 13 + 7F 20 - 6F 21 
+ 6F 22 + 6F 30 - 4F 31 + F 40 ) - 15{N + F 02 - 2F n + F 20 )(N + F 04 + 6F 03 
+ 7F 02 - 2F n - 6F 12 - 4F 13 + 7F 20 - 6F 21 + 6F 22 + 6F 30 - 4F 31 + F 40 ) 

- 6X1 (K x - F 05 - 10F 04 - 25F 03 - 15F 02 + 15F 12 + 20F 13 + 5F 14 + 15F 20 

- 15F 21 - 10F 23 + 25F 30 - 20F 31 + 10F 32 + 10F 40 - 5F 41 + F 50 ). (A6) 
In the above equations all cumulants K n can be directly measured due to the equality: 

F ik = -rr/ifc, (A7) 
PiP 2 

and 

N=(N 1 ) + (N 2 )= { -^+ { -^, (A8) 

Pi P2 

K X = (N\) - (N 2 ) = M _ W . (A9 ) 

Pi P2 
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